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In this work wc use the hcavy-quark-light-diquark picture to study the semileptonic decay 
Ah ^ Ac + 1 + 9i in the so-called hybrid scheme. Namely, we apply the heavy quark effective 
theory (HQET) for larger (corresponding to small recoil), which is the invariant mass 
square of I + v, whereas the perturbative QCD approach for smaller to calculate the form 
factors. The turning point where we require the form factors derived in the two approaches 
to be connected, is chosen near pcut = 1.1. It is noted that the kinematic parameter p 
which is usually adopted in the perturbative QCD approach, is in fact exactly the same as 
the recoil factor uj = v ■ v' used in HQET where v, v' are the four velocities of Ab and Ac 
respectively. We find that the final result is not much sensitive to the choice, so that it is 
relatively reliable. Moreover, we apply a proper numerical program within a small range 
around pcut to make the connection sufficiently smooth and we parameterize the form factor 
by fitting the curve gained in the hybrid scheme. The expression and involved parameters 
can be compared with the ones gained by fitting the experimental data. In this scheme the 
end-point singularities do not appear at all. The calculated value is satisfactorily consistent 
with the data which is recently measured by the DELPHI collaboration within two standard 
deviations. 

PACS numbers: 14.20. Mr, 14.20.Lq, 12.39.Hg, 12.38.Bx, 12.38.-t 



I. INTRODUCTION 

The general theory of QCD has been developed for more than 40 years, and at present, nobody 
ever doubts its validity. However, on the other side there is still not a rehable way to deal with 
the long-distance effects of QCD which are responsible for the quark confinement and hadronic 
transition matrix elements, because their evaluations cannot be done in perturbative approach. 
Thus, one needs to factorize the perturbative sub-processes and the non-perturbative parts which 
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correspond to different energy scales. The perturbative parts are in principle, calculable to any 
order within the framework of quantum field theory, whereas the non-perturbative part must be 
evaluated by either fitting data while its universality is assumed, or invoking concrete models. 

The perturbative QCD method (PQCD) has been applied to study processes where transitions 
from heavy mesons or baryons to light hadrons are concerned , namely the PQCD which 

includes the Sudakov resummation, is proved to be successful for handling processes with small 
4- momentum transfer q^. Indeed the processes involving heavy hadrons may provide us with an 
opportunity to study strong interaction, because compared to Aqcd there exist natural energy 
scales (heavy quark masses) which can be used to factorize the perturbative contributions from the 
non-perturbative effects. On the other hand, for the processes involving heavy hadrons, at small 
recoil region, where v ■ v' is close to unity {v and v' denote the four-velocities of the initial and 
final hadrons), i.e. the momentum transfer is sufficiently large, the heavy quark effective theory 
(HQET) works well due to an extra symmetry SUf{2) ® SUs{2) [4]. Therefore the HQET and 
PQCD seem to apply at different regions of q'^. For a two-body decay, the momentum transfer is 
fixed by the kinematics, however, for a three-body decay, q'^ would span the two different regions. 

Among all the processes, the semi-leptonic decay of hadrons plays an important role for probing 
the underlying principles and employed models because this process is relatively simple and less 
dependent on the non-perturbative QCD effects. Namely leptons do not participate in strong in- 
teraction, and there is no contamination from the crossed gluon-exchanges between quarks residing 
in different hadrons which are produced in the weak transitions, whereas such effects are important 
for the non-leptonic decays. Thus one might gain more model-independent information, such as 
extraction of the Cabibbo-Kobayashi-Maskawa matrix elements from data. In the semi-leptonic 
decays of heavy hadrons it is expected to factorize the perturbative and non-perturbative parts 
more naturally. Recently the DELPHI collaboration reported their measurement on the A;, decay 
form factor in the semi-leptonic process A;, — > A^-l- / + and determined the parameter in the 
Isgur-Wise function ^{v ■ v') = 1 — p"^ {1 — v ■ v') [5|. 

There is a flood of papers to discuss the semi-leptonic decays of heavy mesons and the concerned 
factorization. By contraries, the studies on heavy baryons are much behind P, 3) because 
baryons consist of three constituent quarks and their inner structures are much more complicated 
than mesons. In this work, we are going to employ the one-heavy-quark-one-light-diquark picture 
for the heavy A^ and Ac to evaluate the form factors of this semi-leptonic transition A^ — > Ac + l + v. 
Even though the subject of diquark is still in dispute, it is commonly believed that the quark- 
diquark picture may be a plausible description of baryons ol, especially for the heavy baryons 
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which possess one or two heavy quarks. 

The kinematic region for the semileptonic decay A;, ^ Ac + / + can be characterized by the 
quantity p, which is defined as 

where p, p' are the four-momenta of A5 and Ac respectively. It is noted that this parameter 
p which is commonly adopted in the PQCD approach is exactly the same as the recoil factor 
uj = V • v' used in the HQET. The momentum transfer q'^ in the process is within the range of 
(mi + my)'^ < <f_< (-^Afc — -^Ac)^5 equivalently, it is 1 < p < j^'' = Pmnx ■ In the framework 



of the HQET [J], this process was investigated by some authors [6|, |7|, [lOj]. For larger q'^, the 
HQET works well, whereas one can expect that for smaller q'^, the PQCD approach applies. In this 
work, following Ref.Q], we calculate the contribution from the region with small g^, i.e. p — > Pmax, 
to the amplitude in PQCD. One believes that the PQCD makes better sense in this region. Korner 
et al. discussed similar cases and suggested that the symmetry for smaller recoil is different from 
that for larger recoil, so they used the Isgur-Wise function to obtain the amplitude in the kinematic 
region of smaller recoil, but Brodsky-Lepage function for larger recoils Our strategy is similar 



that we apply the PQCD for small q^ while apply HQET for large where the PQCD is no longer 
reliable, instead. 

Concretely, when integrating the amplitude square from minimum to maximum of q^ to gain the 
decay rate, we divide the whole kinematic region into two parts, small and large {p, equivalently). 
We phenomenologically adopt a turning point at a certain p value, to derive the form factors 
(defined below in the text) in terms of PQCD in the region from pmax to this point and then 
beyond it we use the HQET instead [s]. We let the two parts connect at the turning point. From 
Ref. P], we notice that as p < 1.1, the PQCD result is not reliable, so that we choose the turning 
point at vicinity of p = 1.1. To testify if the choice is reasonable, we slightly vary the values of the 
turning point as choosing p = 1.05, 1.10 and 1.15 to see how sensitive the result is to the choice. 
Moreover, it is noted that as pcut = 1-1 and 1.15 are chosen, the two parts connect almost smoothly. 
Even though, to make more sense, we adopt a proper numerical program to make the connection 
sufficiently smooth, namely we let not only the two parts connect, but also the derivatives from 
two sides are exactly equal. In fact, small differences in the derivatives are easily smeared out by 
the program. Later in the text, we will explicitly show that the final result is not much sensitive 
to it, thus one can trust its validity. We name the scheme as the "hybrid" approach. We also 
parameterize the form factor with respect to p based on our numerical results. In fact, when we 
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integrate over the whole kinematic range of p, we just use the parameterized expression. 

Moreover, we not only reevaluate the form factors fi and gi of the exclusive process in the 

diquark picture, but also calculate the form factors /2 and 52 which were neglected in previous works 
3]. So far, the data on the A;, semi-leptonic decay are only provided by the DELPHI collaboration 
5|] and not rich enough to single out the contributions from /2 and (72 • More accurate measurement 

in the future may offer information about them. Our treatment has another advantage. In the 

pure PQCD approach, there is an end-point divergence at p — > 1, even though it is mild and the 

decay rate which includes an integration over the phase space of final states, i.e. over p, is finite. 

As calculating the contribution from the region with large (p — > 1)^ to the form factors in terms 

of the HQET, the end-point divergence does not exist at all. 

We organize our paper as follows, in section II, we derive the factorization formula for A;, AcW. 

Our numerical results are presented in Section III. Finally, Section IV is devoted to some discussions 

and our conclusion. 



II. FORMULATIONS 



The amplitude of A;, Ac^z/decay process is written as: 



M = -j=Vcbh''{l - 75)z^i < Ac{p ) I C7^(l - 75)^ I Mp) >, 



(2) 



where p and p' are the momenta of A;, and Ac respectively. According to its Lorentz structure, the 
hadronic transition matrix element can be parameterized as 



Mf, = < Ac(p') I C7^(l - 75)6 I Afe(p) > 



Up') 



/2((?') +7552(g') + T7^(/3(g') +7593(9')) 



Ma, y 'I Ma, 

where q=p-p' and cj^,^ = [7^, 7^]/2. 

For the convenience of comparing with the works in literature, we rewrite the above equation 
in the following form according to Ref.js] 



(3) 

Up), 



= Up) 



7m (FM^)+l^G^{q^)) + -^(F2{q')+T.G2{q')) 



+ ^{F3{q')+J,Gs{q')) 

Ma, 



Up)- 



(4) 



^ If nil is not zero, p cannot be exactly 1, thus the superficial singularity does not exist at all, but the form factors 
are obviously proportional to 1/m; which has the singular property. 
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For the case of massless leptons, 

- l^>i = 0, (5) 

thus the form factors fz and result in null contributions. The contributions from /2 and §2 were 
neglected in previous literature Q], nevertheless in our work, we will consider their contributions 
to the matrix elements and calculate them in terms of the diquark picture and our hybrid scheme. 
The kinematic variables are defined as follows. In the rest frame of A;, 

P={P ,P ^Pt)=[-^,^,Ot), (6) 

and 



the diquark momenta inside A;, and Ac are parameterized respectively as 



ki = [0, -^xukirj , k2 = [-^^1^2, 0, k2Tj , (8) 
where p = jz^-jj—,^±^ P~^\/ P^ ~ 1) ^^i = /p^ and X2 = fcj/p'^. According to the factorization 



theorem 



la ], the hadronic matrix element is factorized in the b-space as 

= dxidX2 J j^i^j^;^'ifA,[X2,h2,p 

xHf,{xi,X2,hi,h2,r, Ma^ , *At (xi , bi , j5, /i) , (9) 



where r = Ma^/Ma^. The renormalization group evolution of the hard amplitude is shown as 
follows Q 

P 



i?^(xi,X2,bi,b2,r,MA,,/i) = exp[-4 / — 7g(a5(;u))]iJ^(a::i, 2:2, bi, b2, r, Ma^, t<j(fe)),(10) 

J u P 



where 7g(as(/x)) is the anomalous dimension. 

The wave function of A5 which has the heavy-quark and light-diquark structure, is given as 

*A,(xi,bi,p,/i) = /f^$f^(xi,bi,^j,/x)xA,Afe(p,A), (11) 

where Afe(j>, A) is the baryon spinor, and the superscript S denotes scalar diquark (spin= 0, isospin= 
0). jf^^ is a constant introduced in literature. Xa^ flavor component of the baryon, namely 

XAj — ^'^^\u where 6^ and S]^^ ^ are the creation operator of b-quark and the scalar diquark 

of quarks. 
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FIG. 1: The lowest order diagrams of the hard parts of the transition processes in the quark-diquark picture. 

The Ac distribution amplitude b similar form, 

*A,(2;2,b2,p',/i) = /£$f^(x2,b2,p',Ai)xiAc(p',A2). (12) 

Including the Sudakov evolution of hadronic wave functions, i.e. running the scale of wave function 
from n down to u)i{ui2) [l^ : 

A^{xi,hi, p, n) = exp[-s{u}i,{l - xi)p~) - 2 / — 75(0^(71) )]^>A(,(xi, bi), 
$i^(x2,b2,y,/x) = exp[-s(u;2,(l-X2y+)-2 T ?79(«s(7^))]^i(^2, ba), (13) 

JuJ2 

where uJi = l/&j(i=l,2). 

In our work, the effective gluon-diquark vertices are defined as [ijj, ll 

SgS: igse{pi+p2)^Fs{Q^), 

2\ _ X ^0 



S.=a,{Q^)la.(Ql) (if Q^>Ql): S. = l {if Q^<Ql), (14) 

where = —{pi — P'lf' ■ 

According to the factorization scheme which is depicted in Figdl it is straightforward to obtain 
the analytic expressions of the form factors Fi, Gi, F2, G2, -F3 and G3, by comparing eqs. ([4]) with 
([9]). In the above derivation, the following transformation has been used. 

kx=Ap^Bpl, k2 = Cp + Dp (15) 

and the explicit expressions of A,B,C and D are given in the appendix. 

Then we can obtain the analytical form of fi and gi (i=l,2,3) making use of the relations between 
them and Fj, Gj(j=l,2,3) listed below 

1 F2 F3 2 1. G2 G3 x,^2 

I F2 ^3 -.,^2 1/^2 Gs 2 f-ia\ 
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We do not display the expressions of /a and 53 for the reason given above. The form factors are 
integrations which convolute over three parts [l^: the hard-part kernel function, the Sudakov 
factor and the wave functions of the concerned hadrons as 

rl roo /'27r 

fiiOi) = ATTCpfKjk, dxidx2 bidbib2db2 dO ^A^{x2,h2) (17) 
Jo JO Jo 

xkernSgSf^(^g^)^Atixi,hi) exp [-S{xi,X2, bi, b2, r, Ma J], 

where the explicit expression of the kernel functions kernSgS f.(^g.-^ is given in the appendix for 
concision of the text. The explicit form of the Sudakov factor appearing in the above equations is 
given in Ref.[l^ as 



5(xi,X2,bi,b2,r,MA^ 

= s{uJi, (1 - xi)p~) + s(a;2, (1 - X2)p'+) + 2 f^^'^ '^^^[as{Ti)) + 2 T"''' ^^^(0,(71)), (18) 



where Cp = 4/3 is the color factor. 
The wave function ^A^{xi,k.iT) is 



'^Nbxl{l-xi)^ 

<fA,{x,MT) - ^[(1 xi)2:e2(1 - x,) + £^,(1 - xrfxj + kj^]^ ' ^^^^ 



and 



^>A,(xi,6i) = y'd2kiT$A,(xi,kiT)e"^^^-^^ 



Nbx\{l - x^YbiK2{J{l -ab- xi)2xf(l - xi) + Epbil - xifxjbi) 



(20) 



2[(1 -ab- xi^xlil - xi) + Epbil - xO^xf] 

where K2 is the modified Bessel function of the second kind. If neglecting the transverse momentum 
kir, i.e. set ki^ ~ 0, the wave function can be simplified as 

U,) ~ Nbxl{l-x^) 

The normalization conditions are set as [3] 

/ <^A^{xi)dxi = 1, 

JO 

/■I A 

jl^..{x,)xldx, = (^ + ^)- (22) 

The first formula determines the normalization of the parton distribution of the baryon, whereas the 
second one is related to the effective mass of the light diquark A ~ M^^, —f^b^ and the third formula 
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reflects connection between hadronic matrix element of the kinematic operator Ai = — 2M/^ I 
by{iDj_)'^by I Ah) and hadronic distribution amplitude. To satisfy the above three normalization 
conditions, the parameters would take the following values A = 0.848 GeV, Ai = —0.76 GeV^, 
Uh = 0.916, Epb = 0.0051 and Nf^ = 0.0219, and in the following numerical evaluation, we will use 
them as inputs. 

For the wave function of Ac, the expressions are the same as that for Af,, while the corresponding 
parameters are A = 0.8849 GeV, Ai = -1.87GeV2 , Uc = 1.48, £pc = 0.080 and = 12.14. Thus 
we can write the differential decay width as 
dV _ Gl\V,b\^ 



1 



with 



dp 247r3 

[ - l)(3r2 - Apr + 2r + 3) + \gi\'^{p + l){3r'^ - Apr - 2r + 3) 

+6 /i /2(r +l){p- l)(r2 - 2 rp + 1) + 6 51 52(r - l){r^ -2rp+l) 
+ \f2?{p - l)(3r2 - 2pr + 4r + 3)(r2 - 2pr + 1) 

+ |52p(/0 + l)(3r2 - 2pr - 4r + 3)(r2 _ 2pr + 1) ], (23) 



l<P<Prna. = \r4P^+'4P^]. (24) 



2 \ Ma, Ma, 



III. NUMERICAL RESULTS 



A. The results of PQCD 



In the one-heavy-quark-one-light-diquark picture, the {ud)^ diquark in A;, is considered as a 
scalar of color anti-triplet. To calculate the form factors in the framework of PQCD, one can 
adjust the product /a^/ac *° empirical formula fi{p 

max) ~ l'32/yO^(j2, given by the authors 
of Ref.[3|]2. Then the corresponding parameters are obtained as 

rrik ~ Ma, = 5.624 GeV, ~ Ma, = 2.2849 GeV, V^b = 0.040, 

/aJa. = 0.0096GeV2, ^ 3.22GeV2, Aqcd = 0.2GeV. (25) 

Using these values, we can continue to numerically estimate the form factors /i and gi. Figl2^a) 
shows that the form factor /i is exactly equal to \gi\ in the heavy quark limit. The form factor /2 



When first calculating the form factors, there were no data available, the authors of [3 used a reasonable estimate 
of BR{Ab — > Acliy) as about 2%. Nowadays, measurements have been done and with the data, we have made a 
new fit. 



9 



1.8-1 




0.0 -L, , 1 , 1 I I I . 1 -T ' 1 ' 1 ' 

1.0 1.1 1.2 1.3 1,4 1.5 ^ 

(a) ^ (b) 



FIG. 2: (a) Form factors /2, /i and \gi\ (b) Differential decay width of A;, — > AJiy 

and g2 are much smaller than /i and l^il, thus they can in fact be safely neglected for the present 
experimental accuracy. 

In Figl2jb), we plot the dependence of the differential width dV /dp on p. Although both /i and 
gi have an end-point divergence at p ^ 1 in PQCD approach, the differential decay rate is finite. 

If one extrapolates the PQCD calculation to the region with smaller p values, we obtain the 
form factors /i^2 and gi^2 in that region where there are obvious end-point singularities at /) — > 1. 
Redo the computations with the extrapolation (in the original work Q|, the authors extend the 
tangent of the PQCD result at a small p value to p = 1, so the end-point singularity is avoided) 
and obtain BR(Ah — > AJi') which is about 1.35% (slightly smaller than the value of 2% guessed 
by the authors of Ref.^], because then no data were available. ). 

Obviously, the calculation in PQCD depends on the factor /a^/a^, which regularly must be 
obtained by fitting the data of semileptonic decays, so that the theoretical predictions are less 
meaningful. Instead, we will use our hybrid scheme where we do not need to obtain the factor 
/yVft/Ac fitting data, since the connection requirement substitutes the fitting procedure (see below 
for details). 
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B. The results of HQET 

The transition rate was evaluated in terms of HQET by the authors of Ref. [s]. According to 
the definitions given in eq. (fT6|) , we re-calculate /i , /2 , gi, §2 while dropping out /a and and also 
obtain similar conclusion that fi and gi are the same in amplitude, but opposite in sign, as shown 
in Fig. [2]^a), whereas /2 is very small and g2 is exactly zero in HQET. The theoretical prediction 
on the rate of the semi-leptonic decay — > Ac + / + in the HQET is 



T{Ah — > AJi') = 1.54 X 10 GeV, without contributions from /2 and g2, 

r(Ab AcW) = 1.56 X 10^^"^ GeV, with contributions from /2 and g2, (26) 

and the branching ratio is 2.87%. 

The transition rate of A^ AJiy has recently been measured by the DELPHI Collaboration Q], 
and the value of BR{Ai, — > AJi') is 5.0'tl'l{stat)~^\[2{syst)%. It is noted that the result calculated 
in terms of HQET is only consistent with data within two standard deviations. 

C. The hybrid scheme: ReconciUng the two approaches 

As widely discussed in literature, in the region with large (small p— values), the result of 
HQET is reliable, whereas for the region with small (larger p— values) the PQCD is believed to 
work well. Therefore, to reconcile the two approaches which work in different p regions, 1 < p < 
Pmax, we apply the HQET for small p, but use PQCD for larger p. Our strategy is that we let 
the form factors /i, gi derived in the PQCD approach be equal to the value obtained in terms of 
HQET at the point pcut- The numerical results of the form factor fi = \gi\ in the hybrid scheme 
are shown in Fig l3] for three different pcMt— values: 1.05, 1.10 and 1.15, respectively. It is noted 
that for Pcut = 1-1 and 1.15, the left- and right-derivatives are very close and the connection is 
smooth, whereas, as pcut is chosen as 1.05, a difference between the derivatives at the two sides 
of Pcut = 1-05 obviously manifests. We then adopt a proper numerical program to smoothen the 
curve, namely let the derivatives of the two sides meet each other for any p— value near the cut 
point. Figl3] shows that such treatment in fact does not change the general form of the curve, but 
makes it sufficiently smooth for all p values including the selected cut point pcut- 

Another advantage of adopting such a "hybrid" scheme is that there does not exist end-point 
singularity for the form factors at p = 1. Since at the turning point, we let the form factors derived 
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FIG. 3: Form factors /i = in the hybrid scheme 



in terms of PQCD be connected with that obtained in HQET, the product /a,,/Ac is automatically 
determined by the connection. With the value, we calculate the form factors within the range of 
small in PQCD. In this hybrid scheme, one does not need to invoke the data on the semi-leptonic 
decay to fix the parameter /a^/Ac at all. 

By our numerical results obtained in the hybrid scheme, the form factors f\ (or g\) can be 
parameterized in a satisfactory expression, here we only present the expression for p^ut = 1-10 as 

h{p) = 1 - 3.61(/5 - 1) + 7.24(p - 1)2 - 5.83(p - 1)^, (27) 



and similar parameterized form factors were discussed in ref.[17|. 

The expression can be described by only one "Isger-Wise function" for the transition A;, A+ll 
at the heavy quark limit, and it is done by the DELPHI collaboration based on their data on 
Af, — > AcW- It is parameterized as Q] 

^{u;) = l-p\L0-l) + O{{u;-lf), (28) 

where p"^ = 2.03 ± 0A6{stat)tiJ)l{syst). 

Obviously, this expression is only valid to the leading order, i.e. linearly proportional to a; — 1 
where lv exactly corresponds to the parameter p which is commonly adopted in the PQCD language. 
By contrast, our result includes higher power terms because the 1/M corrections are automatically 
taken into account in our work. It is noted that the coefficient of the linear term in our numerical 
result is reasonably consistent with the p^ obtained by fitting data. 

To obtain the total decay width, we integrate over the whole range p from 1 to Pmax, the 
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FIG. 4: The Gi/Fi ratio obtained in our hybrid scheme 



integrand is the parameterized form factor in eq. (l27p . We obtain 



r(A;, ^ AJu) = 1.65 X lO^^^GeV, BR = 3 



with p,ut = 1.10, /aJa, = 0.0149 GeV' 



(29) 



For a comparison, we present the results corresponding to other two pcut values where the smoothing 
treatment is employed, and they are 

r(Ab ^ AJu) = 1.54 X 10^^"^ GeV, BR = 2.87% pcut = 1-05, /a Ja, = 0.0142 GeV^; 
r{Ab AJu) = 1.67 X 10"^^ GeV, BR = 3.12% p^ut = 1-15; /a Ja, = 0.0150 GeV^ (30) 

One can notice that the factor /a^/Ac does not change much as pcut varies and they are about 1.5 
times larger than the value obtained in pure PQCD (eq. (j25p ). When the turning point is chosen at 
Pcut = 1-05, the branching ratio calculated in the hybrid scheme is very close to the result obtained 
in pure HQET, while for pcut = 1.10 and pcut = 1-15, the resultant branching ratio is slightly larger 
than that obtained in pure HQET, but more coincides with the data. 

In our scenario, the HQET is applied for smaller p, and the values of the form factors at pcut 
are fixed by the theory. The values can also determine /a^/Ac which will be used for the PQCD 
calculations for larger p. The HQET is an ideal theoretical framework, but there is an unknown 
function which is fully governed by the non-perturbative QCD effects, that is the famous Isgur-Wise 
function. The function can be either obtained by fitting data, or evaluated by concrete models. 
Various models would result in different slopes. The authors of ref. ^ used the Drell-Yan type 
overlap integrals to obtain the slope which is what we employed to get the parametrization eq. (l27p 
and the slope is —3.6. Instead, the authors of evaluated the slope in the Isgur-Wise function by 
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means of the QCD sum rules. According to their result, we re-parameterize the form factor fi{p) 
and have 

flip) = 1.01 - 1.57(/3 - 1) - 2.59(/> - 1)^ + 6.99(p - 1)^ (31) 
Correspondingly, we obtain 

r(A6 ^ Ac/z^) = 2.38 X lO^i^GeV, Si? = 4.45% with pc«t = 1-10, /a Ja, = 0.0213 GeV^. 

The fitted slope by the DELPHI collaboration is = 2.03 ± 0A6{stat)tiJil{syst) [sl], which 
is between the two theoretical evaluated values. In these references, only linear term remained, 
due to uncertainties in the approximations the deviations are understandable. Therefore, one 
can note that there is a model-dependence which mainly manifests in the slope of the Isgur-Wise 
function. Even though they deviate from each other at the linear term, the high power terms 
would compensate the deviation slightly and the predicted values on the branching ratio in two 
approaches eqs. (|27l3ip are qualitatively consistent with data. There indeed is a byproduct which 
brings in an advantage that more accurate measurements can help to make judgement on validity 
of the models by which the Isgur-Wise functions are evaluated. 

To make more sense, we purposely present the ratio Gi/Fi obtained in the hybrid scheme in 



Fig. U it is noted that the ratio is qualitativel y c onsistent with that given in ref. ll| which was 



shown on the left part of Fig. 3 of their paper 
The authors of ref. 
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Ul] extended p into the un-physical region [p > p max for Afe — > Kflu), while 
we only keep it within the physical region. I t is noted that in the physical region, numerically our 
result is very close to that obtained in ref. 



lll |. But if one extends the curve to larger /?, he will 
notice that our curve is convex, but theirs is concave, namely the coefficient of the quadratic term 
has an opposite sign, but the difference is too tiny to be observed or bring up substantial difference 
for the evaluation of the decay width. 



IV. DISCUSSION AND CONCLUSION 



In this work, we investigate the semi-leptonic decay ^ Kc + I + I'l in the so called "hybrid 
scheme" and the diquark picture for heavy baryons A;, and Ac. The hybrid scheme means that 
for the range of smaller (larger />— value) we use the PQCD approach, whereas the HQET for 
larger (near p = 1), to calculate the form factors. We find that the form factors /2 and g2 
can be safely neglected as suggested in the literature. Besides, we do not need to determine the 
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phenomenological parameters /a^ and /a^ by fitting data in tlie hybrid scheme as one did with the 
pure PQCD approach. Our result is generally consistent with the newly measured branching ratio 
Afe — > Ac+l+i>i within two standard deviations and the end-point singularities existing in the PQCD 
approach are completely avoided. In fact, the final result somehow depends on the slope in the 
Isgur-Wise function of the HQET, which is obtained by model-dependent theoretical calculations. 
Therefore, we may only trust the obtained value to this accuracy, the further experimental data 
and development in theoretical framework will help to improve the accuracy of the theoretical 
predictions. 

The quark-diquark picture seems to work well for dealing with the semi-leptonic decays of 
A;,, and we may expect that the quark-diquark picture indeed reflects the physical reality and is 
applicable to the processes where baryons are involved, at least for the heavy baryons 18|]. This 



picture will be further tested in the non-leptonic decays of heavy baryons. We will employ the 
diquark picture and PQCD to further study the non-leptonic decay modes in our future work. 
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APPENDIX A: EXPLICIT EXPRESSIONS OF HARD KERNEL 



kernSgSf, = { [{BC - AD){1 + r) - {C^ - A{2 - C + Dr) + {B + D)r{Dr - 2p) + 



kernSgS, 



91 



kernSgSf^ 



C{Br + 2Brp + 2Drp - 2)) 



{BC - AD)r (1 + r) - {A^ + A 



hha + 

(C - r{D + 2p- 2Bp)) + r{B'^r - 2{Dr + Cp) - B{C + (-2 + D)r 
+2Cp)))]hhb}MlFs{r^), 

{ [{BC - AD){r -1) + (C^ + A{-2 + C + Dr) + {B + D)r{Dr - 2p) 



C{rB - 2Drp - 2rBp + 2)) 



{BC - AD) r(l -r) + {A'^ + A 



hha + 

{C + r{D -2p + 2Bp)) + r{B'^r - 2{Dr + Cp) + B {{-2 + D)r - C 

+2Cp)))]hhb] MlFs{r]), 

{AD - BC) {hha + r hhb) MlFs{r]), 



kernSgSg^ = {AD - BC) {hha - r hhb) mI^Fs{v), 



with 



hha = as{ta)Ko{JxiX2CiMA^MA^bi)Ko{Jx2CiMA,MA^\hi+h2\), 



ta = max{Jx2CiMA^MA,, l/|bi + bsj, 1/&1, 1/&2) 



tb = maxU xiCiMa^Ma,, l/|bi + b2|, 1/62)- 
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(Al) 



where 

77 = Ma^Ma^xiXs^i. (A2) 

The expUcit expressions of A, B, C and D 

^ ^ xi 6 ^ 2:1 Ma, 



6-6 Ma,(6-6) 

JWA6(42-4lj ?i-?2 



The expUcit expressions of hha, hhb are 



hhb = asitb)Ko{^xiX2^iMA,MAM)Ko{\/xiCiMA^MA,\hi + h2\), (A5) 

with 
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